Ordered magnetic phases of the frustrated spin-dimer compound Ba 3 Mn 2 8 
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Ba3Mn2C>8 is a spin-dimer compound based on pairs of S = 1, 3d 2 , Mn 5+ ions arranged on a triangular 
lattice. Antiferromagnetic intradimer exchange leads to a singlet ground state in zero-field. Here we present the 
first results of thermodynamic measurements for single crystals probing the high-field ordered states of this ma- 
terial associated with closing the spin gap to the excited triplet states. Specific heat, magnetocaloric effect, and 
torque magnetometry measurements were performed in magnetic fields up to 32 T and temperatures down to 20 
mK. For fields above H c i ~ 8.7 T, these measurements reveal a single magnetic phase for H ||c, but two distinct 
phases (approximately symmetric about the center of the phase diagram) for H _L c. Analysis of the simplest 
possible spin Hamiltonian describing this system yields candidates for these ordered states corresponding to a 
simple spiral structure for H\\c, and to two distinct modulated phases for H _L c. Both single-ion anisotropy 
and geometric frustration play crucial roles in defining the phase diagram. 

PACS numbers: 75.45.+j, 75.40.-s, 75.30.-m, 75.50.-y 



I. INTRODUCTION 

Antiferromagnetic exchange on a triangular lattice leads to 
geometric frustration - the system cannot satisfy all of the pair- 
wise interactions simultaneously, such that the minimum en- 
ergy of the system does not correspond to the minimum en- 
ergy of all local interactions. The classical solution to the 
Heisenberg antiferromagnet on a triangular lattice with only 
nearest neighbor interactions is the well known 120° struc- 
ture. In this case, the main effect of the frustration is simply 
to produce a non-colinear structure. In an early attempt to 
find the groundstate of the two-dimensional quantum triangu- 
lar antiferromagnet, Anderson proposed the Resonating Va- 
lence Bond (RVB) spin liquid state 1 -.. Subsequent theoretical 
work has indicated that for the simple case with only near- 
est neighbor interactions the classical solution is in fact sta- 
ble against quantum fluctuations, but with a much reduced or- 
dered moment 2 ^. Other models corresponding to more com- 
plex lattices and interactions are still the subject of intense 
theoretical investigation 4 * In this paper, we experimentally 
examine the slightly more complex case of a triangular lattice 
decorated by vertical magnetic (spin) dimers, realized by the 
compound BasMi^Os. 

Spin dimer compounds comprise pairs of strongly coupled 
magnetic ions. Antiferromagnetic intradimer exchange leads 
to a ground state that is a product of singlets, but an applied 
magnetic field can be used to close the spin gap to excited 
triplet states, resulting at low temperatures in a state charac- 
terized by long range magnetic order (LRMO)^. Close to the 
critical field, the effective Hamiltonian that describes the low- 
energy degrees of freedom of such a system can be expressed 
in terms of an effective spin- \ model, or equivalently a lattice 
gas model of hard-core bosons -. The exchange anisotropy in 
this effective model, which is sensitive to the lattice geometry 
and contains contributions from both the interdimer coupling 
and also other anisotropies present in the system, is equiva- 
lent to the balance between potential and kinetic energies in 
the bosonic picture. This anisotropy plays a critical role in 



determining the nature of the ordered stated. 

Spin dimer compounds offer several specific advantages 
over simple (non-dimerized) magnetic lattices. First, varia- 
tion of an external magnetic field can be used to tune the triplet 
density, providing easy access to more of the quantum phase 
diagram, and of course the quantum critical point (QCP). The 
effects of quantum fluctuations at the QCP can be dramatic, 
especially for spins arranged on a frustrated lattice. For ex- 
ample, in the case of the spin dimer compound BaCuSi2C<6, 
the effects of "order from disorder" are suppressed at the QCP 
precisely because the size of the moment is tuned by the exter- 
nal field (i.e. this is an amplitude-driven QCP), and the frustra- 
tion implicit in the body-centered tetragonal lattice of this ma- 
terial ultimately leads to a form of dimensional reduction^' 9 . 
Second, the interdimer exchange energy "protects" the system 
against symmetric anisotropies, effectively suppressing the in- 
fluence of interactions such as dipolar coupling and, for the 
case in which the kinetic energy of the bosons dominates, en- 
abling a realization of a Bose Einstein condensate (BEC) in a 
temperature range which is even comparable to the anisotropy 
energy 10 . (This is why one specifically looks for such an effect 
in a spin dimer compound and not a simple antiferromagnet.) 
And third, spin dimer compounds provide a means to "engi- 
neer" large exchange anisotropies, providing access to lattice 
gas models that would otherwise require unphysical parame- 
ters for a simple antiferromagnet. Specifically, although in- 
dividual exchange couplings in such materials may be nearly 
isotropic ( J z ~ J X y), in the strong coupling limit (intradimer 
exchange J 3> interdimer exchange J') the effective Hamil- 
tonian derived from perturbation theory and acting on the sin- 
glet and triplet states can have strongly anisotropic effective 
exchange couplings (J%*' ^ '"). As an example, in some 
circumstances a large uniaxial anisotropy (JJ 3> J^y) 
would provide favorable conditions for realization of a spin 
supersolid phase-! 1 -. 

In this paper we examine the high field behavior of the 
novel 5 = 1 triangular spin dimer compound BaaM^Os- 
The interest in this specific compound stems from the pos- 
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sibility to explore the interplay between geometric frustration 
(which typically favors uniform spiral states), and single ion 
anisotropy (favoring specific orientations of the moments) in 
the context of a spin dimer compound, in which case we can 
explore the entire quantum phase diagram, ultimately includ- 
ing the quantum critical behavior. We determine the magnetic 
phase diagram via heat capacity, magnetocaloric effect and 
torque magnetization measurements of single crystals, find- 
ing a single magnetic phase for fields aligned parallel to the 
crystalline c-axis, but two magnetic phases for fields oriented 
away from this direction. Comparison of these results with 
a low energy effective model indicate that the competing ef- 
fects of interdimer coupling and single ion anisotropy on the 
triangular lattice stabilize complex magnetic structures char- 
acterized by multiple independent order parameters, including 
spontaneous bond order. 

Ba3Mn2C»8 crystallizes in the rhombohedral R3m structure, 
and is comprised of pairs of Mn 5+ 3d 2 S = 1 ions arranged 
vertically on hexagonal layers (see fig. []I a ))^- Each Mn ion 
is coordinated by distorted oxygen tetrahedra, which results 
in an orbitally non-degenerate ground state. Successive layers 
are stacked following an 'ABC sequence, such that the dimer 
units on adjacent planes are positioned in the center of the tri- 
angular plaquets of the layers above and below. Pairs of ions 
on each dimer are coupled antiferromagnetically, resulting in 
a singlet ground state (fig. [TJb)). Initial powder inelastic neu- 
tron scattering (INS) indicated an intradimer exchange energy 
Jo ~ 1.61(3) meV-Si The same INS measurements revealed a 
spin gap of 1 .05 meV and a hierarchy of additional exchange 
energies, in which interdimer coupling between Mn ions in the 
same plane is characterized by J2 — J3 = 0.112^qqq3 meV 
and interdimer coupling between Mn ions residing on adjacent 
planes is characterized by J\ = — 0. 062^066 me V- Prelim- 
inary single crystal INS measurements further refined these 
values to Jo = 1.65meV, J2 — J3 = 0.109meV and Ji = 
— 0.120meV— (Here we preserve the labelling of exchange 
energies initially suggested by Uchida et al^-). BagM^Os 
can therefore be described as a quasi-2D material in which 
planes of vertical dimers arranged on triangular layers inter- 
act weakly in the perpendicular direction^. Electron para- 
magnetic resonance (EPR) experiments in the diluted com- 
pound Ba3(Vi_ a; Mn 2 ;)208 (where the V 5+ 3d ion carries no 
moment) reveal a single ion uniaxial anisotropy characterized 
by D = 5.81 GHz, corresponding to 0.024 meV— . Similar 
measurements for the pure compound BaaM^Os indicate a 
zero field splitting of the triplet characterized by \D\ = 0.032 
meV— , the difference being due to the presence of additional 
symmetric anisotropics in the concentrated lattice, in particu- 
lar dipolar coupling between the two ions on each dimer. Val- 
ues of the g-tensor were revealed by EPR measurements of the 
diluted compound to be g cc = 1.96 and g aa = g^ = 1.97—, 
and are confirmed by EPR measurements in the concentrated 
lattice^. 

The spin gap in BaaM^Os can be closed by an applied 
field H c i ~ 9 T, and measurements of powder samples have 
revealed a magnetization which rises approximately linearly 
with field from H c i until the eventual triplet saturation field 
is reached at H C 2 ~ 26 T— . Heat capacity measurements of 
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FIG. 1 : (Color online) (a) Schematic diagram showing arrangement 
of MnC>4 tetrahedra in R3m structure of BaaM^Og. Intradimer (Jo), 
out-of-plane interdimer (Ji), in-plane direct interdimer (J2) and in- 
plane crossed interdimer (J3) exchange bonds are drawn as thick 
black, dotted black, thin black and dashed black lines, respectively^. 
Only two of the in-plane crossed interdimer exchange bonds (J3) 
were drawn for clarity, (b) Energy spectrum as a function of field for 
an isolated dimer composed of 5=1 Mn 5+ ions with antiferromag- 
netic exchange J. (c) Susceptibility of BaaMifeOg as a function of 
temperature for fioH = 0. IT applied parallel (up triangle) and per- 
pendicular (down triangle) to the c-axis. Mol refers to one dimer 
unit. Red line shows fit to dimer model as described in main text. 

polycrystalline samples revealed tantalizing evidence for two 
phase transitions 19 , but to date no single crystal samples have 
been available to definitively determine the phase diagram, 
nor probe the ordered states. A second spin gap to excited 
quintuplet states can be closed by increasing the field beyond 
H C 2, leading to an additional increase in magnetization above 
approximately 30 T— . In this paper we investigate solely the 
nature of the long range magnetic order which results from 
closing the spin gap to the triplet states. 

II. EXPERIMENTAL METHODS 

Single crystals of BaaMr^Os were grown by a slow cool- 
ing flux method. Polycrystalline precursor was synthesized 
by mixing MnO and BaCOs reagents and sintering at temper- 
atures up to 1000° C for 100 hours with intermediate regrind- 
ings according to the previously published method 13 . To grow 
crystals of BaaMr^Os from solution requires a strongly oxi- 
dizing flux - we have found anhydrous NaOH to work espe- 
cially well. Polycrystalline material was mixed with NaOH in 
a molar ratio of 1 to 19 and placed in a 20 cc alumina crucible 
lightly sealed with a cap. The mixture was heated to 700° C 
over the course of 24 hours, allowed to dwell for 5 hours, and 
then slowly cooled to 300° C at which temperature the furnace 
was turned off. The flux can be removed by repeated washes 
in water. Crystals grown by this method form as hexagonal 
tablets, and have a mass of up to 100 mg. From initial X-ray 
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diffraction measurements they have a small mosaic spread of 
less than 0. 10°, and minimal magnetic impurities (see charac- 
terization in following section). 

Low-field susceptibility measurements were performed us- 
ing a commercial Quantum Design MPMS XL SQUID mag- 
netometer for fields of 1000 Oe applied both parallel and per- 
pendicular to the c-axis. 

Heat capacity (C p ) data were collected down to 0.35 K us- 
ing a thermal relaxation time technique, both in a Quantum 
Design Physical Properties Measurement System (PPMS) for 
fields up to 14 T applied parallel to the a and c axes, and at 
the National High Magnetic Field Laboratory (NHMFL) us- 
ing a home-built calorimeter in a resistive magnet for fields up 
to 32 T applied parallel to the a-axis. Results were repeated 
for several crystals. Magnetocaloric effect (MCE) scans were 
also performed in the resistive magnet at the NHMFL using 
the same home-built calorimeter used for heat capacity mea- 
surements. Data were taken for sweep rates of 2 T/min and 5 
T/min, for both increasing and decreasing fields, and for fields 
oriented both parallel and perpendicular to the c-axis. 

Cantilever torque measurements were performed at the 
NHMFL for fields up to 1 8 T using a superconducting magnet. 
The sample was mounted on a cantilever forming one plate 
of a capacitor and oriented such that the magnetic field was 
slightly misaligned with one of the principle axes. As a con- 
sequence of the misalignment and the intrinsic g-anisotropy, 
the resulting magnetization was not parallel to the applied 
field, resulting in a finite torque. Although the g-anisotropy 
of B^M^Os is relatively weafeii, the resulting torque is 
nevertheless measurable and can be used to provide a sensi- 
tive probe of changes in the magnetization close to H&, and 
in particular at the phase transition(s). The cantilever was 
mounted on the cold finger of a dilution refrigerator, and field 
sweeps were performed for temperatures from 20 to 820 mK. 



III. RESULTS 

The low field susceptibility of B^M^Os is shown in fig. 
Q2c) for fields aligned parallel and perpendicular to the c-axis. 
There is negligible anisotropy, consistent with the small g- 
anisotropy determined from EPR measurement s 17-18 . As pre- 
viously described for polycrystalline samples 13 , the temper- 
ature dependence of the susceptibility can be well fit to an 
isolated 5=1 dimer model: 



(a) H 1| c 



Xd = 



2N{3g 2 f i 2 B (1 + 5e- 2 ^ J ) 
3 + eP J + 5e- 2 ^ J 



(1) 



if one includes a meanfield correction to account for ex- 
change with neighboring dimers of the form \ = i+\ Xd 
where A = 3 [Ji + 2 ( J 2 + J 3 )] / {Ng 2 fi 2 B ) and (3 = l/k B T. 
The fit also includes terms to account for a temperature- 
independent background (xo) and a small concentration of 
paramagnetic impurities (C/T). These fits result in estimates 
of J = 1.44±0.01 meV and g c ~ g a ~ 2.01±0.03, which 
are close to the values obtained from INS 14 and EPR' 8 experi- 
ments, and also to similar susceptibility fits for polycrystalline 
samples' 3 . The fit is rather insensitive to the precise value of 
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FIG. 2: (Color online) Heat capacity (shown as C p /T) as a function 
of temperature for fields applied (a) parallel and (b) perpendicular to 
the c-axis. Blue (red) lines indicate the transition into Phase I (Phase 
II) from higher temperatures as determined by peaks in C V IT as a 
function of field. 



A (which is why this is a poor method to estimate interdimer 
exchange coefficients), but nevertheless returns a best value 
of A = 5.0± 0.3 mol/emu which is remarkably close to the 
calculated value of 6.6 mol/emu based on estimates of the ex- 
change constants 3 J\ + 6( J2 + J3) = 0.83 meV obtained from 
single crystal INS measurements^ 4 -. The T-independent term 
has a value xo = 2xl0~ 4 emu/mol, and the impurity Curie term 
corresponds to just 0.4% unpaired Mn 5+ ions. 

Heat capacity data are shown in fig. [2] for temperatures 
down to 0.35 K. For fields greater than H c \ applied parallel 
to the c-axis (fig. |2ja)) there is just one phase transition in 
this temperature range. We label this ordered state Phase I 
to distinguish from a distinctly different phase observed for 
fields parallel to the a-axis. The data show a lambda-like tran- 
sition, suggestive of 3DXY ordering, with an integrated en- 
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tropy which increases as the field is increased from H c \ up to 
the maximum field for which data were taken for this orienta- 
tion, which was 14 T. These data points are included as solid 
square symbols in the phase diagram shown in fig. [3 a). 

For fields oriented parallel to the a-axis, heat capacity data 
were taken up to H C 2 ~ 26 T (fig. |2b)). These data show a 
remarkable sequence of phase transitions at low temperature 
with an unusual division of entropy. For fields between 9 and 
1 1 T just one transition is observed above 0.35 K; between 1 1 
and 13 T two distinct transitions are clearly resolved; for in- 
termediate fields, only a single transition; between 24 and 25 
T two transitions are again observed; and finally close to the 
triplet saturation field only one transition is observed. Anoma- 
lies in the heat capacity marking these phase transitions are 
joined by red and blue lines in fig. |3b), and T c values in- 
cluded in the phase diagram shown in fig. [^b) as solid sym- 
bols. The heat capacity anomaly for phase transitions joined 
by the blue line in fig. |2jb) are lambda-like, similar to those 
observed for fields parallel to the c-axis, and accordingly we 
label this state Phase I. The integrated entropy associated with 
these transitions first rises with field, and then after H is in- 
creased beyond the midpoint of the phase diagram, reduces in 
magnitude again. In contrast, the anomaly associated with the 
phase transitions connected by the red lines in fig. |2jb) are 
less divergent, and although the data do not permit a critical 
scaling analysis, nevertheless are more suggestive of an Ising 
transition. We refer to this state as Phase II. The rise in C p /T 
associated with this phase transition does not appear to vary 
with field within the uncertainty, indicating that the change in 
entropy is only weakly dependent on the applied field. Esti- 
mates of the integrated entropy are difficult due to the close 
proximity of the second phase transition, and also due to the 
large background magnetic contribution to the heat capacity 
associated with the other gapped states. However, a crude es- 
timate of this entropy was obtained for several fields for which 
the only resolvable transition is between the disordered phase 
and Phase II (plotted in fig. He) for H = 10.5 T). Upper 
and lower bounds for the entropy were determined by assum- 
ing a minimum and maximum possible background, shown as 
blue and red lines respectively in fig. |4|c), yielding an average 
of 0.45 ± 0.20 J/molK. Within the uncertainty, this value ap- 
pears to be symmetric for fields above and below the midpoint 
of H c i and H C 2 as shown in fig. |4jd), and for this reason we 
also label the ordered state on the right hand side of Phase I in 
fig. Etb) as Phase II. 

Additional insight to these phase transitions is provided by 
MCE measurements. Phase transitions are evident from a 
sharp increase (decrease) in the temperature of the sample on 
entering (leaving) the ordered state. In practice, points on the 
left (right) hand side of the phase diagram (open symbols in 
fig. [3]) were determined from a sharp peak (trough) in the 
first derivative of the temperature with respect to field taken 
on up (down) field sweeps, each corresponding to the case of 
entering the ordered state. These data are in close agreement 
with heat capacity measurements, with small differences be- 
ing ascribed to differences in sample alignment and, where 
two different calorimeters were used, thermometry. A slight 
assymmetry in the magnitude of the change in temperature 




FIG. 3: (Color online) Phase diagram obtained from heat capacity 
(solid squares) and MCE (open triangles) measurements for fields 
applied (a) parallel and (b) perpendicular to the c-axis. MCE traces 
are shown in gray (green) for decreasing fields for sweep rates of 2 
(5) T/min. Representative data for increasing fields are designated by 
arrows and shown in blue for sweep rates of 5 T/min in panel (a) and 
2 T/min in panel (b). Labels indicate Phases (I) and (II), as described 
in the main text. Dashed line in panel (a) indicates the anticipated 
phase boundary for fields between 16 and 21 T based on similar data 
in panel (b). Shaded vertical bands are guides to the eye to draw 
attention to the broad features observed in MCE measurements cen- 
tered at 8.8T and 26.5T. 



between up and down sweeps provides evidence that phase 
boundaries between the disordered state and Phase I, between 
the disordered state and Phase II, and between ordered Phases 
I and II are all weakly first order in this temperature range. 

The total change in entropy associated with a phase transi- 
tion can be calculated from MCE measurements by summing 
the increase in entropy of the system plus the entropy lost as 
heat to the bath from the sample stage: 



6S t = - 



C (T i+1 - Tj) + k {Tj - T base ) 
Ti 



(2) 



where i labels successive temperature points taken as the field 
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FIG. 4: (Color online) Entropy estimates for the Phase II for fields 
oriented parallel to the a-axis. (a) Three representative MCE mea- 
surements (dashed lines, right axis) and the corresponding change in 
entropy, SSi, (solid lines with the same color, left axis), (b) Total en- 
tropy associated with phase transitions seen in MCE measurements, 
calculated as described in the main text. Horizontal arrows indicate 
the entropy associated with the lower transition, (c) Upper and lower 
bounds of the integrated entropy associated with the phase transition 
observed in heat capacity at 10. 5T. (d) Entropy on entering Phase II 
from MCE (open circles) and C p (solid squares). 



is swept (typically 40 evenly spaced points per Tesla) and k is 
the thermal conductivity of the thermal link in the calorime- 
ter. A linear interpolation for k was calculated as a function 
of temperature and field for the calorimeter, and values of 
the heat capacity were taken from measurements performed 
in the PPMS calorimeter. Figure |4ja) shows three represen- 
tative MCE data sets for fields oriented perpendicular to the 
c-axis (dashed line, right axis), and the associated change in 
entropy SSi between successive data points (solid curves, left 
axis). All three data sets were taken for increasing fields and 
for the same sweep rate of 2 T/min. As can be seen, SSi 
shows two successive peaks as a function of field, which cor- 
respond to the two phase transitions. A practical measure 
of the change in entropy associated with each phase transi- 
tion is therefore provided by the integrated entropy up to the 
minimum in SSi, which is shown in fig. SJb). For the low- 
est temperature data set (black curves) the two transitions are 
well separated and the total entropy associated with the first 
transition exhibits a clear plateau. For the higher temperature 
sweeps (red and green curves) the two transitions are slightly 
closer in field, and the total entropy exhibits more of a kink 
than a plateau. Nevertheless, these data allow an estimate of 
the integrated entropy associated with each transition, which 
are plotted in fig. |4jd) for these and some additional interme- 
diate temperature sweeps. Within the uncertainty of this anal- 
ysis, the change in entropy associated with entering Phase II 
from the disordered state is essentially independent of temper- 
ature, with an absolute value that agrees remarkably well with 
the value extracted from heat capacity measurements (square 
symbol in fig. SJd)). In contrast, the change in entropy asso- 



FIG. 5: (Color online) Field and temperature dependence of the mag- 
netic torque (scaled by the field strength) for fields applied close to 
perpendicular to the c axis. A linear interpolation scheme between 
field sweeps taken at 20, 45, 75, 114, 160, 183, 237, 280, 318, 372, 
412, 417, 470, 535, 600, 696 and 800 mK (dark lines) has been 
used to generate the 3D surface. Red and blue lines indicate phase 
transitions determined as described in the main text. Labels indicate 
Phases I and II. 



ciated with entering Phase I depends strongly on temperature, 
consistent with inspection of the heat capacity data shown in 

fig. Kb). 

The MCE experiments did not reveal any evidence for ad- 
ditional phase transitions other than those indicated in fig. [3] 
However, broad features were observed centered at H c \ and 
H C 2 for both field orientations (vertical shading in fig. [3]). This 
effect is associated with the rapid change in magnetization 
with temperature at these fields. Specifically, Maxwell's rela- 
tion ^jr\ H P = ^jj\t,p implies a large change in the entropy 
at H c i and the triplet saturation field for field sweeps per- 
formed at constant temperature, or conversely, a large change 
in the temperature of the sample in a constant entropy ex- 
periment. The effect reverses sign for increasing/decreasing 
fields, and for entering/leaving the field region between H c \ 
and the triplet saturation field, and becomes rapidly smeared 
out at higher temperature. 

The phase diagram obtained from heat capacity and MCE 
measurements (fig. |3j is almost but not completely symmetric 
about its midpoint. For fields oriented parallel to the a-axis 
(fig- EJb)), Phase II is slightly narrower in field on the right 
hand side of the phase diagram relative to the left. The ab- 
sence of strict particle-hole symmetry is consistent with the 
large bandwidth relative to the spin gap observed in powder 
INSH. 

Further information about the magnetically ordered phases 
can be gained from torque magnetometry. Raw data taken for 
fields nearly perpendicular to the c-axis, scaled by the mag- 
netic field strength, are shown as a 3D surface in fig. The 
sharp increase in torque at 8.5 T for the lowest temperatures 
corresponds to the spin gap closing as the minimum of the 
the S z =l triplet band crosses the singlet leading to a finite 
magnetization. As temperature is increased, thermal effects 
smear this rapid increase in torque, and field derivatives of 
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FIG. 6: Representative cantilever torque measurements for fields ap- 
plied close to perpendicular to the c-axis, expressed as t/H and its 
first two derivatives for temperatures of 535 mK (a-c) and 75 mK 
(d-f). Vertical dashed lines mark phase transitions as described in 
text 
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t/H at H c i rapidly broaden with temperature (figures [6] and 
[7]). Superimposed on top of the broad rise in torque with field, 
the two phase transitions seen in heat capacity and MCE ex- 
periments are clearly visible as breaks in the slope of t/H, 
and points on the phase boundary can be extracted from field 
derivatives as described below. 

Representative torque data and their derivatives are shown 
in fig. [6ja-c) and (d-f) for temperatures of 535 and 75 mK re- 
spectively. Specifically, the data in panels (a-c) are represen- 
tative of field scans for temperatures above 160 mK, whereas 
panels (d-f) are representative of field scans below 160 mK. 
Considering first the 535 mK data, both phase transitions are 
characterized by a decrease in the slope of t/H above each 
critical field (fig. |6ja)). The first derivative of t/H (fig. |6tb)) 
more clearly reveals this effect, exhibiting a relatively broad 
minimum close to each transition. Inspection of fig. |6ja) also 
reveals a distinct downwards "hook" in t/H at both critical 
fields, being more prominant for the second transition. Ex- 
trapolation of t/H from above and below the two transitions 
reveals a discontinuous (downwards) change at the critical 
field, indictive of a first order transition. In this case, an em- 
pirical estimate of the critical field can be found in the sharp 
inverse peak in the second derivative (fig. |6|c)) which effec- 
tively marks the onset of the change in magnetization (dashed 
line in fig. [SJ. From simple energetic reasons the magneti- 
zation should always rise as the field is increased, therefore 
these features indicate a change in the anisotropy of the or- 
dered states at the critical fields. The absence of a delta-like 
transition in the heat capacity indicates that both transitions 



FIG. 7: (Color online) (a) First and (b) second derivatives with re- 
spect to field of torque divided by field for fields close to perpendic- 
ular to the c axis, (c) Phase diagram extracted from MCE (triangles), 
C p (squares), and torque (circles). Open symbols are used to signify 
a weakly 1 st order transition. Closed symbols signify a second or- 
der transition. Red (blue) symbols indicate the lower (upper) field 
transition. 



are only weakly first order. 

Below 160 mK, typified by data shown in figures HJd-f), 
the upper transition is still marked by a distinct decrease in 
the slope of t/H as the field is increased, and features in the 
two field derivatives at this transition are similar to those seen 
at higher temperatures in fig.s[6ja-c). The data do not exhibit 
as clear of a "hook", but the smooth evolution of the raw data 
and both derivatives from the behavior at high temperatures 
implies that this transition is still weakly first order. How- 
ever, for temperatures below 160 mK the lower field transition 
is qualitatively different to the upper transition. Specifically, 
the torque exhibits a rapid increase in slope (panel (d)), lead- 
ing to an abrupt step in the first derivative (panel (e)) and a 
sharp peak in the second derivative (panel(f)). This behavior 
is consistent with a second order phase transition, in which 
case consideration of the critical exponents associated with 
this thermal phase transition implies that the sharp positive 
peak in the second derivative marks the critical field (vertical 
dashed line)2i. 

Figures[7](a) and (b) show 3D surface plots of first and sec- 
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ond field derivatives of t/H respectively for all temperatures 
measured, demonstrating the smooth evolution of the various 
features described above. Points on the resulting phase di- 
agram, extracted as described above, are shown in fig. |7Jc) 
and are in agreement with points taken from MCE and C p 
measurements up to the inherent angular misaligment in this 
torque measurement. For temperatures near 160mK (i.e. in 
the temperature range at which the phase transition changes 
from 2 nd to 1 st order) there is some ambiguity as to which 
feature marks the lower field transition, and in these cases 
values obtained from both criteria are plotted. Unfortunately 
there are insufficient data points below 160 mK to extract a 
meaningful critical scaling exponent for the phase boundary 
approaching the QCP. 

Careful inspection of fig. [5]reveals a small rise in the torque 
at approximately 6 T for the lowest temperatures. No sharp 
features are associated with this rise, indicating the absence 
of any ordering transitions, and it is rapidly smeared out with 
increasing temperature (fig. |7ja)). The observation of a fi- 
nite magnetization for fields below U c \ indicates the presence 
of terms in the spin Hamiltonian that mix singlet and triplet 
states, such as antisymmetric Dzyaloshinskii-Moriya (DM) 
interactions. The midpoint of the dimer unit in BaaM^Os 
is a center of inversion symmetry, ruling out the presence of 
an intra-dimer DM term. However, DM interactions are still 
possible on the inter-dimer bonds, both within planes and be- 
tween adjacent planes. The absence of any torque at low fields 
indicates that such terms are relatively weak. Further EPR 
measurements should clarify the origin of this feature. 



IV. DISCUSSION 

In the absence of measurements that directly determine the 
magnetic structure of the ordered phases of BasM^Os, we 
resort to an analysis of the spin Hamiltonian which describes 
the system: 



H 



j2 ^s^-s^+D^Ky 



ifi 



(3) 



Here i, j designate the coordinates of the dimers while \i,v = 
1, 2 denote each of the two spins on a given dimer. The intra- 
dimer exchange interaction is Jq = Jm2 = Ji2i\- The nearest- 
neighbor (NN) inter-dimer exchange interaction on the same 
layer is J2 = an< ^ ^ = w i tn I 1 ^ v > where the 

dimers i, j belong to the same layer. Finally, the NN interlayer 
interaction is given by Ji = J^ji, where i, j denote the po- 
sition of NN dimers on adjacent layers. The quantization axis 
(z-axis) is always along the field direction. The components 
a = a,b,c and r\ = x, y, z depend on the field orientation; i.e. 
a = c, 77 = z for H\\c, and a = a, b, r\ = x for H _L c. We 
use values for the exchange parameters determined by inelas- 
tic neutron scattering 1 *'! 5 , an( j values for D and the g-tensor 
determined by EP R 17 ' 18 , as described in the Introduction. The 



midpoint of the magnetic dimer is a center of inversion sym- 
metry, so we need not consider an intradimer Dzyaloshinskii- 
Moriya (DM) interaction. The effect of intradimer dipolar 
coupling is implicitly contained within the anisotropy term 
D (S v ) 2 together with the single ion anisotropy. For this ini- 
tial analysis we assume that interdimer dipolar and DM in- 
teractions are negligible, although some subtle features of the 
torque magnetization may require inclusion of these terms for 
a more complete description. 

Since Jq ^> J±, J2, J3, we can include the inter-dimer 
terms perturbatively relative to the intradimer Heisenberg 
term. For this purpose, we keep only the low energy singlet, 
1 00), and triplet |11) states of the single dimer problem for 
H ~ H c i . We use a pseudospin s = i to represent the low- 
energy singlet and triplet states: |00) — and |11)— *|f). The 
low-energy effective Hamiltonian to first order in J\ , J2 , J3 
and D, H, results from projecting H into the low-energy sub- 
space generated by the states |00) and |11): 
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where each dimer has been replaced by an effective site, 
B = g aa n B H - J - 3 ( J 2 + J 3 ) /2 - 3Ji/4 - 
I is the layer index, indicates that the sites i and j are 
nearest neighbors (NN) on the same layer and denotes 
NN on adjacent layers. In addition we have a(z) = (the 
model is U(l) invariant for H\\c) and a(x) = —8D/3Jq (see 
Appendix A). Although the exchange anisotropy terms are of 
higher order ( 0(JiD / Jq) with i = 1,2) than the rest of the 
terms, we need to include them to give account of the observed 
differences between the H \\ c and H _L c cases. 

To understand this effective Hamiltonian, and some proper- 
ties of its ground states, it is simplest to consider the two ex- 
perimentally observed cases - fields applied parallel and per- 
pendicular to the c-axis. In each case we first consider the 2D 
lattice, effectively setting J\ = (i.e. no interplane coupling), 
and then we consider the full 3D case. Finally we compare 
how the interlayer frustration is relieved in BaaM^Os while 
the interlayer frustration is preserved in the similar 2D spin 
dimer compound BaCuSi2C»s. 



A. Fields parallel to c 

Considering first the 2D lattice (effectively setting J\ = 
0), and noting that the D anisotropy does not act on the s xy 



components of the pseudospins (i.e. a(z) = 0), the system 
consists of independent triangular layers of vertical dimers, 
and the effective Hamiltonian reduces to 



(J2-J3) 
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(g ccf i B H - Jo - 3 (J 2 + J 3 ) /2) £ sf, (5) 



The effective exchange anisotropy is easy-plane, i.e., the XY 
component of the exchange dominates. At T = 0, the triplets 
condense (canted XY antiferromagnetic ordering) for H > 
H c i = (Jo — 4 ( J 2 — J3)) / (<7 cc /xb) into a state that can be 
approximated by a direct product of single dimer states of the 
form: 
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The canting angle 
netic field, while $ 
The relative phase between different layers is determined by 
4>l, which can take any value for the moment because we are 
assuming that Ji = 0. The expectation values of the pseu- 
dospins take a simple form 



= 9 is uniform and is set by the mag- 
= <t>i + Q • Ti with Q = (±± ±\)a*. 
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This corresponds to a canted antiferromagnetic state, in which 
the s xy component of the pseudospins orient 120° with re- 
spect to each other to minimize the interdimer exchange en- 
ergy (i.e. the s xy components of the pseudospins on each tri- 
angular plaquet sum to zero) analogous to the classical solu- 
tion for a Heisenberg AF on a triangular lattice. This is illus- 
trated in fig. [HJa). The phase, corresponding to the angle of the 
s xy component of the pseudospins relative to the crystal lat- 
tice, spontaneously breaks the U(l) symmetry of the effective 
Hamiltonian, and the ordered state can be described as a Bose- 
Einstein condensate. In terms of the original spins on each Mn 
site, the ordered state still consists of a canted AF with the s xy 
component of the spins on adjacent dimers oriented 120° with 
respect to each other, but with these components reversed for 
spins on the top and bottom of each dimer unit. 

Now consider the 3D lattice with nonzero J\. The Hamil- 
tonian still maintains U(l) symmetry, and the ordered state 
will still correspond to a triplet condensate (consistent with the 
lambda anomaly observed in heat capacity measurements for 
this field orientation) because the XY interaction in the Hamil- 
tonian dominates the Ising interaction. However, the system 
now has the possibility to gain additional energy from the in- 
terlayer exchange. The s xy component of the total spin on any 
triangular plaquette for the classical case in fig. Oa) is zero, 
but if the pseudospins twist around the z-axis to form a spiral 
structure in which successive spins along the [100] and [010] 
directions rotate by an angle a = 120° ± e in the XY plane, 




[010] 



[110] 



[100] 



(b) a<120° 
q=0 




(c)a>120° 
q z =n J ; <0 




FIG. 8: (Color online) Schematic diagrams of the predicted spin 
structure for field applied along c-axis. (a) 120° structure on a trian- 
gular plaquette corresponding to the classical solution for Heisenberg 
spins on a 2D lattice. Green and red arrows indicate s 2 and s xy com- 
ponents of the pseudospin representing each dimer unit, respectively. 
Inset shows crystal axes. Full 3D structure of s xy components of 
pseudospins on successive layers for values of a less than and more 
than 120° leading to ordering wavevectors along the z-direction of 
and ic* are plotted in (b) and (c) respectively. Dashed blue arrows 
indicate the total pseudospin moment on each triangular plaquette, 
illustrating ferromagnetic interplane coupling. 



the system is then able to benefit from the interlayer coupling. 
There are two degenerate solutions that minimize the total en- 
ergy. The first solution is characterized by a uniform phase 
along the c-axis: tfii = 0. In contrast, the phase is staggered 
4>i = lir, in the second solution. In addition, the shift of the 
single-layer ordering wave-vector from Q = (±i, ±-|)a* to 
Q = ±(a, a )^f (note that a* = -^=- in this non-orthogonal 
basis) has opposite signs for the two cases: 



cos a 



■/1 



4(J 2 - J 3 ) 
Ji 

4(J 2 - J 3 ) 



for 



= 



for cf>i = In, 



(8) 



corresponding to angles a = 111° and a = 129° respec- 
tively for the values of Ji and ( J 2 — J3) obtained from single 
crystal INS. The second solution was previously reported in 
Ref. [HI The gain in energy of this in-plane twisting due to 
the inter-layer interaction, J\, scales linearly with e (i.e. for a 
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given phase relation between adjacent layers there is a "right 
way" and a "wrong way" to twist the spiral structure). In con- 
trast, the loss in intralayer energy from breaking the perfect 
120° structure scales quadratically with e because a = 120° 
is the minimum energy structure for J\ = 0. Hence, an arbi- 
trarily small interlayer exchange J\ is able to stabilize a spiral 
phase with an incommensurate wave vector as indicated by eq. 
(|8). The resulting phase is incommensurate, which has been 
confirmed by preliminary single crystal nuclear magnetic res- 
onance measurements and neutron scattering measurements at 
high fields^. 

A similar structure, but with antiferromagnetic interplane 
coupling, has already been proposed by Uchida et al. fol- 
lowing their initial estimation of the exchange constants in 
BasM^Ogii However, two subtleties to the ordered phase 
were not anticipated in that earlier paper. The first of these is 
that e can take both positive and negative values because there 
are two degenerate solutions: <fii = or <fii = In. If a < 120°, 
as illustrated in fig. [8jb), then for J\ < the component of 
the ordering wavevector along the c direction, q z , is equal to 
(fa = 0). However, if a > 120°, as illustrated in fig. [SJc), 
then the component of the ordering wavevector along the c di- 
rection, q z , is equal to ^c*) (<fri = In), leading to a doubling 
of the unit cell along the c-axis. A second subtlety of the 
ordered phase is that the resulting structure contains triangu- 
lar plaquettes on which pair of spins are more closely antifer- 
romagnetically aligned along specific directions, while other 
pairs of spins are less perfectly antiferromagnetically aligned 
along equivalent crystallographic directions. Specifically, ad- 
jacent spins along the [110] direction have a relative angle 
of 2a = 2(120° + e) = 120° - 2e, in contrast to adjacent 
spins along the [100] and [010] directions which have a rela- 
tive angle of a ~ 120° + e. Magnetostriction associated with 
this "bond ordering" must lead to a subtle lattice deforma- 
tion, characterized by a separate but related order parameter. 
If the coupling to the lattice is strong enough, we can even 
anticipate that the phase transition will become weakly first 
order—. Microscopic energetics of the two different bond or- 
dering scenarios will break the degeneracy and determine the 
c-axis ordering wavevector. 



B. Fields perpendicular to c 

For fields oriented away form the c-axis, the anisotropy 
term D breaks the U(l) symmetry of the Hamiltonian, and 
can stabilize an Ising-like modulated structure. To understand 
the nature of this phase, it is instructive to first consider the 
case in which D vanishes on a 2D lattice ( J\ =0) and a mag- 
netic field is applied perpendicular to the crystalline c-axis, for 
instance along the direction midway between [100] and [110] 
(i.e. along the "point" of a triangular plaquette). All of the 
same arguments given above for the case H\\c and J\ = still 
apply, but the quantization axis now lies in the ab plane so the 
X Y order lies in the plane defined by the two vectors [00 1 ] and 
[010] (fig. |9ja)). Since D = 0, there is no anisotropy in this 
plane, and the pseudospins spontaneously break U(l) sym- 
metry - the ordered state is a Bose Einstein condensate. How- 



Hie 

(a)D=J,=0 (v-0) 




(b) D<0 - Phase I (0<Y<n/2) 




(c)D<0 -Phase II (y^tt/2) 




FIG. 9: (Color online) Schematic diagram showing spin structure for 
fields pointing along the 'point' of a triangular plaquette obtained 
by minimizing parameters in eq. [9] Green and red arrows indicate 
s z and s xy components of the pseudospin representing each dimer 
unit, respectively, (a) Spin structure for D=0 and Ji =0 for H \ | [1 10] . 
The s xy components of the pseudospins are oriented 120° from each 
other, equivalent to the case shown in figure[8{a) for H ||c, but rotated 
into the [010]-[001] plane, (b) Partially modulated phase (Phase I) 
for D < and Ji finite. The s xy components of the spin precess 
along an elliptical path with an incommensurate wavevector close to 
(3 , ^)a*. (c) Maximally modulated phase (Phase II) stabilised close 
to H c i and H C 2- 



ever, a finite value of D in equation|4]qualitatively changes the 
nature of the ground state. For the specific example of the field 
oriented along the tip of the triangular plaquette, a negative 
value for D, appropriate for BasM^Os, implies an easy axis 
for the s xy component of the the pseudospins along the [001] 
direction. For the 2D lattice (i.e. J\ — 0), the in-plane interac- 
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tion between spins still favors a magnetic structure for which 
the total spin in the XY plane vanishes. To minimize both 
the anisotropy energy and also the in-plane exchange energy, 
the system may adopt an inhomogeous magnetic structure in 
which the component of the pseudospins along the hard axis 
are depressed relative to along the easy axis while the canting 
angle along the z-axis is adjusted so as to preserve zero net 
s xy spin on each triangular plaquette (fig.HJb)). 

Considering now the 3D lattice (finite Ji), the general form 
of the pseudospins describing this modulated behavior in a 
spin dimer system with uniaxial anisotropy is 



(4>u\s x \4>u) 
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where < 7 < ir/2 sets the ratio between the maximum 
amplitude of the s v and s x components and consequently the 
amplitude of the modulation of the s z component (7 = is 
unmodulated and 7 = ir/2 is maximally modulated). The 
positive (negtive) sign in the last line of eq. [9] holds for 
9 > 7r/4 (9 < 7r/4). Minimization of the Hamiltonian © 
with respect to the different parameters 9, 7, Q and <fii yields 
a ground state. The optimal values of Q and <f>i are still very 
well approximated by eq. dS). The optimal values of 9 and 
7 as a function of the field H are shown in figfiO] The or- 
dered ground state has no s v component of the pseudospin 
(7 = 7r/2) for H near H c \ (fig. |9jc)). However the effec- 
tive exchange anisotropy in H (eq. H} penalizes the modula- 
tion of the s z component and favors a less modulated struc- 
ture (7 < 7r/2) when the z-component of the real magne- 
tization becomes large enough, i.e., when H — H c \ is large 
enough. This is presumably the origin of the two distinct 
phases observed in thermodynamic measurements for fields 
oriented perpendicular to the c-axis, and indeed a full anal- 
ysis including all three triplet states quantitatively accounts 
for H c i—. Even without including these terms the agreement 
with the measured phase diagram is remarkable. According 
to our results shown in fig. [10] the transition between both 
phases is of second order at T = 0. For H slightly larger 
than H1.2 (critical field for the transition between phases I and 
II), we obtain the field dependence ir/2 — 7 oc y/H — 
characteristic of a mean field transition. Correspondingly, the 
total magnetization and 9 exhibit a kink at Hi t 2- The result- 
ing structure stable at higher fields is still modulated along the 
z direction but with a finite component along the y direction 
(fig. |9jb)). The s y component of the pseudospin varies to the 
s x component as a function of field, yielding an unmodulated 
structure exactly at the middle of the dome since s z =0 at this 
field, equivalent to the H\\c structure. At this field, rotation of 
field into the H\\c direction therefore occurs without crossing 
a phase boundary, consistent with our labelling of Phase I in 
fig- Ha) an d (b). The energy associated with the anisotropy 
a(x) (J2 — J3) ~ D (J2 — J3) / Jo ~25mK is small, consis- 
tent with the observation of a A-like anomaly in heat capacity 




10 11 12 13 14 15 16 



FIG. 10: Variational parameters 7 and 9 and normalized magnetiza- 
tion M/M S at that result from minimizing the energy (Ti) for the spin 
configuration described by eq. The vertical dashed lines mark 
H c i and the transition between Phase II and Phase I. We have used 
values of J = 19.2 K, Ji = -1.39 K and J 2 - J 3 = 1.26 K that 
result from preliminary fits of triplet dispersion measured at H — 
in a single crystal of BaaM^Og^. The value of J 3 = 1.27 K was 
chosen to fit the measured optimal field: (H c i + H C 2) /2 ~ 17.3 
T. The disagreement between the calculated H c \ ~ 9.6 T and the 
measured value of ~ 8.7 T is due to our two-level (singlet-triplet) 
approximation. Good quantitative agreement is obtained when the 
other two triplets are included^. 



seen in fig. |2jb) (i.e the critical scaling associated with the 
Ising phase transition will only apparent very close to T c ). 

The resulting modulated structures are characterized by 
three separate order parameters: there is a finite modulation 
of the magnetization along the field direction (Ising order), 
there is also Ising ordering along the easy direction perpen- 
dicular to the field, and finally there is bond order. Referring 
back to the heat capacity data (fig. [2j, the anomaly associated 
with this phase transition into Phase II is clearly different to 
the A anomaly seen for fields oriented parallel to the c-axis, 
although the data do not permit a critical scaling analysis. Es- 
timates of the change in entropy through the transition into 
Phase II vary only weakly with field or temperature, implying 
a significant lattice contribution. 



C. Interplane frustration: Comparison to BaCuSi20s 

It is instructive to compare the proposed magnetic struc- 
tures of Ba3Mn2C>8 with those of the closely related spin 
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® j f =rm 5 (f) j, = -(j 2 -4)/io 
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FIG. 1 1 : (Color online) Comparison of the effect of interlayer cou- 
pling on the dispersion minimum for BaCuSi20e (panels a-c) and 
BaaM^Og (panels d-f). (a) Idealized body centered tetragonal lat- 
tice of spin dimer compound BaCuSi20e showing the intralayer 
exchange J' and the interlayer exchange J/, (b) and (c) show 
BaCuSi2C>6 triplet dispersion minimum at (1/2, 1/2) for interlayer 
exchange J j equal to and 0.1J', respectively, (d) Triangular lat- 
tice of Ba3Mn208 showing intralayer exchanges J2 and interlayer 
exchange J\. J3 not shown for clarity (see fig. |TJa)). (e) BasIVk^Og 
triplet dispersion for the case Ji=0, showing the minimum located at 
(1/3,1/3). (f) Including a small interlayer coupling, shown here for 
the specific case of Ji = — (J2 — J-i)/lQ, shifts the minimimum of 
the dispersion away from (1/3,1/3). 



where we have used the same labeling of the in-plane inter- 
dimer coupling (J 7 ) and the inter-plane coupling (Jf) intro- 
duced in ref&. Significantly, since Jf is an order of mag- 
nitude smaller than J', the presence of finite interlayer ex- 
change does not move the minimum of the dispersion from 
Q = (§,|)a* (fig. [TTTb)). The resulting interlayer frustra- 
tion of the perfect bet lattice is partially lifted by the effect of 
quantum fluctuations^, resulting in a fully three dimensional 
magnetic structure away from the QCP. However, since the 
quantum phase transition is driven by amplitude fluctuations, 
the QCP itself is two-dimensional. Additional terms in the 
Hamiltonian describing the real material (including deviations 
from the perfect bet lattice 27 , and the effects of the higher en- 
ergy triplet states—) will partially lift the perfect frustration 
even at T = 0. Nevertheless, the observation in this material 
of 2D critical scaling exponents close to H c i— indicates that 
there is a considerable region over which the fluctuations are 
determined by the 2D fixed point. 

The model describing Ba^Mi^Os has markedly different 
properties. In this case, the minimum of the dispersion rela- 
tion associated with a single triangular layer, derived from the 
effective Hamiltonian (eqn. @J, is located at Q = (i, -|)a* 
referred to the primitive unit vectors of the lattice'-. An ex- 
pansion of the single particle dispersion about the point Q 
now contains terms linear in k: 

E(Q + k)-E(Q) ~ ( J 2 - J 3 ) (hi + k 2 y + k x k y ) 

/3 

+ ^yJl(k x + k y ). (11) 

The competing effects of the linear and quadratic terms means 
that arbitrarily small interplane coupling J\ shifts the min- 
imum of the dispersion from Q = (-|, -|)a*, illustrated in 
fig- ITTT e.f). This has the dual effect of establishing an incom- 
mensurate spiral structure at low temperatures, and lifting the 
perfect interlayer frustration of the 120° structure. In effect, 
the triangular lattice is less frustrated than the bet lattice, at 
least with regard to the effect of interlayer coupling! Hence, 
for Ba3Mn2C<8 we anticipate a fully three dimensional mag- 
netic structure down to T = 0, and three dimensional criti- 
cal scaling exponents. Experiments to determine these are in 
progress. 



dimer compound BaCuSi20g. Both materials have a small 
frustrated interlayer couplin g 14,26 . In the case of BaCuSi2C<6, 
vertical dimers are arranged on a (nominally^ 7 -) body centered 
tetragonal (bet) lattice 2 *. The dispersion of the triplet associ- 
ated with a single square layer of dimers, derived from the ef- 
fective Hamiltonian for this system (see eq.(2) in ref. |28|) . has 
planar inversion symmetry about the 2D minimum, which is 
located at Q = (§, |)a* (fig. [TTJi))- Therefore, an expansion 
of the single particle dispersion about the point Q including 
the weaker interlayer coupling only contains terms quadratic 
in k: 

E (Q + k) - E (Q) ~ J'(k% + k 2 y ) + 2J f k x k y (10) 



V. CONCLUSION 

In summary, via a combination of heat capacity, MCE and 
cantilever torque magnetometry measurements, we have es- 
tablished the low temperature phase diagram of the 5 = 1 spin 
dimer compound Ba 3 Mn208 which is associated with closing 
of the spin gap to excited triplet states. These data reveal two 
distinctly different ordered states (labelled as Phases I and II) 
for fields oriented perpendicular to the crystalline c-axis. For 
fields oriented parallel to c, only Phase I is observed, at least 
for temperatures above 300mK. Analysis of the data, and con- 
sideration of the minimal spin Hamiltonian that describes the 
system, indicates that Phase I consists of a spiral structure 
which is stabilized by the weak interlayer coupling. For fields 
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perpendicular to c this phase is partially modulated due to the 
single ion anisotropy. In contrast, Phase II appears to be a 
fully modulated structure with no moment along the y direc- 
tion, stabilized by a combination of the single ion anisotropy 
and the interlayer coupling. Both phases implicitly contain 
spontaneous bond ordering, and are characterized by multiple 
independent order parameters. Experiments are in progress 
to directly determine the magnetic structures of which we de- 
scribed the qualitative features. 
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APPENDIX A: DERIVATION OF THE EFFECTIVE 
EXCHANGE ANISOTROPY 

Here we derive the effective exchange anisotropy term of 
H that is induced by the single-ion anisotropy term Hsi = 
D^ li M (S^) 2 of TL. There are two different processes that 
contribute to the amplitude, 2a(x) (J2 — J3), of the effective 
exchange anisotropy. In the first process, Hsi induces an in- 
tradimer transition between the singlet 1 00) and the quintuplet 
1 22) and the corresponding matrix element is: 

(22|/J[(^) 2 + (^) 2 ]|00) = -|. (Al) 

The gap between the singlet and |22), for a field H such that 
the (| 11)) triplet and the singlet states are degenerate, is Jo- 
The second step is a transition between the state 1 22) ^ 1 00) j on 
nearest-neighbor dimers i, j and the state with two 5f + 5*| = 



1 triplets |ll)i|ll)j produced by the interdimer Heisenberg 
interactions J2 and J3. The corresponding matrix element is: 

(11^(111^ +F Ja |22) i |00) j = 2(J ^ J3) (A2) 

There is a factor of 2 that results from the fact that the quin- 
tuplet state can be created in any of the two dimers (i or j) 
involved in this second order process. 

The second contribution to 2a(x) ( J2 — J3) comes from the 
following second order process: the pair of singlet states can 
be excited into a pair of triplets with opposite Sf + Sf , 

(lll^llM^ + Hj 3 |00) l |00) J = ^' 2 3 3; , (A3) 
and the triplet state with Sf + S% = — 1 can be flipped into 
the Sf + SI = 1 by the D term, 

(lT|/J[(Sn 2 + (^) 2 ]|ll) = |. (A4) 

Again there are two of these processes because the Sf + 5| = 
— 1 triplet can be in any of both dimers. 

The sum of these contributions leads to an effective ex- 
change anisotropy: 

2 (J 2 - J 3 ) a(x) £ (44 - 44) < A5 > 

with 

a(x) = ~TT- (A6) 

By following a similar procedure, we obtain an effective ex- 
change anistropy for the inter-layer coupling: 

Jra(x) Wi-*a* V ji)- (A7) 

These are the two anisotropic terms that appear in the effective 
Hamiltonian Ti of eq.©. 
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